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Summary. — The isotropic luminosity function (LF) and formation rate history 
(FRH) of long GRBs is by the first time constrained by using jointly both the 
observed GRB peak-flux and redshift distributions. Our results support an evolving 
LF and a FRH that keeps increasing after z = 2. We discuss some interesting 
implications related to these results. 

PACS 98.70. Rz - gamma ray bursts. 



1. Introduction 

The primordial and most extensive information on GRBs obtained directly from obser- 
vations is the differential peak-flux, P, distribution, the so commonly called logN — logP 
diagram (NPD). The NPD is the convolution of basically three factors: the intrinsic LF, 
the FRH and the cosmic volume. In the past, several attempts were done to constrain 
the LF from fits to the NPD, assuming arbitrary GRB FRHs and cosmologies. The re- 
sults were rather poor because the complicated mixing between the model LF and FRH 
introduces a high degeneracy among these factors in the NPD. The most direct way to 
infer both the LF and the FRH is based on the observational luminosity-redshift diagram 
|101 |BJ ^] . Unfortunately, this strategy is so far limited by the small number of GRBs 
with known or inferred z, by the bias introduced by the given method to infer z, and 
by the sensitivity limit related to the z estimate. We propose @] a new strategy for 
constraining the GRB LF and FRH based on the joint use of the NPD and the observed 
or inferred GRB differential z distribution (the N — z diagram, NZD). 



2. The Method and the Data 

The differential NPD and NZD are modelled by seeding at each z a large number 
of GRBs with a given rate, pgrb, and LF, /(Li SO )dLi SO , and then by propagating the 
flux of each source to z = 0. We use the popular ACDM cosmology with f2 m =0.29, 
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f^A = h = 0.71. Here Lj so in the rest frame is defined as L iso = J 30kcV ° ES(E)dE, where 
S(E) is the Band energy spectrum with the (average) parameters taken from [Sj. 
The break energy at rest, Eb, is assumed either constant (512 keV jjj) or dependent on 
L iso 12 (E b = 15 (L iso /10 50 erg s" 1 ) 5 keV). The sensitivity band at z = is fixed to 
the 50-300 keV range of BATSE, . 

We explore two models for the LF, the single and double power laws (SPL and DPL, 
respectively), and two cases, one with a non evolving LF, and another one where Li so 
scales with z as (l + z) s 0H2|. For Pgrb, we adopt the bi-parametric star formation 
rate (SFR) function, psp(z;a,b), given in [7] multiplied by a normalisation factor K, 
and by a function r)(z; c) that allows to control, through its parameter c, the growth or 
decline of Pgrb at z > 2. The strategy is to constrain seven parameters (3 of LF, 3 of 
FRH and S, for details see 0]) by applying a joint fit of the model predictions to the 
observed NPD and NZD (including their errors) . The fitting is based on an extension of 
the Levenberg-Marquardt method to find the least total \ 2 = Xnp + Xnz ■ 

We use the widest sample to date for the NPD It consists of a collection of (i) 
3255 BATSE GRBs longer than 1 s, and with a P limit appropriately extended down 
to 0.1 photons cm~ 2 s _1 , and (ii) a sample of bright Ulysses GRBs with P up to ~ 300 
photons cm -2 s _1 . For the NZD, we use two samples. One consists of a set of 220 
BATSE GRBs with z's inferred by using the luminosity-variability empirical relation 
0, and the other comprises 33 GRBs with known z. The latter is corrected for several 
selection effects. 

3. Results and Discussion 

We have run several combinations of models, with an SPL or DPL LF, adopting an 
Eb constant or dependent on L; so , and either without evolution (6 = 0) or including 
S among the parameters to optimise. The obtained LFs and FRHs are shown in Fig. 
1 (curves) together to the observational data (error bars). A detailed analysis of the 
fittings of these different models to the NPD and NZD data is presented in 0]. In the 
following, we highlight some of the results from this analysis: 

GRB LF and jet angle distribution. As seen in Fig. 1 models with non-evolving LFs 
(SPL or DPL, dotted curves) give in general poor fits to the NPD and NZD. The best 
joint fits are obtained for models with evolving LFs (oc (l + z) s , solid lines). The optimal 
values we find for S are 1.0 ± 0.2. Note that the increasing of Li so with z increases the 
probability to observe GRBs from very high z's. The fits are slightly better for the 
(evolving) SPL LFs than for the DPL ones. The best range of slopes of the SPL LF is 
7 = 1.57 ± 0.03 (/(Lj so )dLi SO oc Li SO ~ 7 ). If the LF is related to collimation effects, we 
have compared our results with the universal structured jet and with the quasi-universal 
Gaussian structured jet models 0]. Our results imply an intermediate case between 
the universal structured jet model with e{9) oc 6~ 2 and the quasi-universal Gaussian 
structured jet model. For the uniform jet model, the jet angle distribution covers an 
indicative range between 2° and 15° at z = 1. These results correspond to the case of a 
luminosity-depending Eb in the model LFs, but we find that for Eb =const the differences 
are not significant. 

GRB FRH and the connection with the cosmic SFR. The models that best fit the 
observed NPD and NZD imply not only an evolving LF but also FRHs, pgrb(-z), which 
steeply increase (by a factor of ~ 30) from z = to z w 2 and then continue increasing 
gently up to z ~ 10 as (1 + z) 1A and (1 + z) for the SPL and DPL LFs, respectively. This 
behaviour of the GRB FRH is qualitatively similar to the one of the cosmic SFRH. Fig. 
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Fig. 1. - Top panel: Peak flux differential distribution (NP) with the axis in the top-right 
part, and z differential distribution (NZ) with the axis in the bottom-left part, both for a SPL 
LF. Error bars show the NP data from and the NZ data according to 5], respectively. 

Dotted lines are for models without evolution, while solid lines are for models with the evolution 
parameter 8 optimised. Thin and thick lines identify the cases with Eb=511 keV or Eb depending 
on Li so , respectively. Dashed straight line is the -3/2 uniform distribution (Euclidean) behaviour 
in the NPD. Bottom panel: Same as in top panel but for a DPL LF. 



2 shows a comparison of the cosmic SFRH traced by rest-frame UV luminosity (corrected 
by dust obscuration) with the GRB FRH translated to SFRH under the assumption of 
a non evolving Initial Mass Function (IMF). The shaded area corresponds to the range 
of GRB models (including the uncertainties) with evolving (8 = 1) SPL and DPL LFs 
that best fit the NPD and NZD data. 

From z = 2- 3toz = 6, the cosmic SFRH traced by UV luminosity decreases 
slightly or strongly 0, while the SFRH linked to GRBs keeps increasing. This dif- 
ference, if proved in future studies, could be related to two effects, (i) A significant 
contribution to the cosmic SFRH at z > 2 seems to come from sources emitting strongly 
in the rest-frame far infrared/submillimetre (not seen in UV). GRBs could be also tracers 
of the dust-enshrouded SFRH of these galaxies (e.g., |||). (ii) The IMF in low-metallicity 
gas (high z's) could be biased toward higher masses compared to the present-day IMF, 
favouring in this way an increasing GRB formation rate with z. From z — 1 — 2toz — 
the UV-luminosity SFRH decreases by less than an order of magnitude, while the SFRH 
linked to GRBs decreases roughly by a factor of 30. Due to the uncertainty on the z 
determination any definitive conclusion about this difference has to be supported by a 
more extended sample with know z for NZD. 

Implications for the progenitors. Our best models give a true (after collimation effect 
correction) GRB FR of ~ 5 10 yr _1 for the Milky Way. Based on astronomical 
arguments we argue that such a FR is close to that of close binary systems consisting of 
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Fig. 2. - Comparison between the observed cosmic SFRH traced by the rest-frame UV luminosity 
from (dots with error bars) and from 2 (triangle) and the SFRH obtained from GRB FRH 
properly normalised (shaded area, see text) 



a WR star and a possible massive BH, with periods of hours. These systems are able to 
generate a massive Kerr BH after the SNIb/c explosion of the WR (helium) star. The 
observational counterparts of these potential GRB progenitors should be luminous X-ray 
binaries (e.g., Cyg X-3), which are estimated to be only a few at present in the MW. 

Implications for cosmology. The finding of the LF properties will shed new light on 
the connection between the true (collimation corrected) energetics of GRBs and their 
spectral features. This connection, studied for the first time by Ghirlanda et al. |13II14| . 
makes GRBs powerful standard candles which hopefully will allow to explore in few years 
the geometry and the kinematics of the universe beyond z = 10 |151 1161 IT7| . 
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